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Inspired by the discovery of the dibaryon d* and the experimental search of NQ dibaryon with 
the STAR data, we study the strange dibaryon NQ further in the framework of quark delocalization 
color screening model and chiral quark model. We have shown NQ is a narrow resonance in AH 
D-wave scattering before. However, the A-H scattering data analysis is quite complicated. Here we 
calculate the low-energy NO. scattering phase shifts, scattering length, effective range and binding 
energy to provide another approach of STAR data analysis. Our results show there exists an ATI 
’’bound” state, which can be observed by the N-Q correlation analysis with RHIC and LHC data, or 
by the new developed automatic scanning system at J-PARC. Besides, we also find that the hidden 
color channel-coupling is important for the NQ system to develop intermediate-range attraction. 
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I. INTRODUCTION 

The possible existence of dibaryon states was first pro¬ 
posed by F. J. Dyson and N. Xuong [l] in 1964. However, 
this topic got considerable attention only after R. Jaffe’s 
prediction of the H particle in 1977 [1]. All quark models, 
lattice QCD calculations and other methods, predict that 
in addition to qq mesons and q 3 baryons, there should be 
multiquark systems: ( qq ) 2 , q^q, q 6 , quark-gluon hybrids 
qqg , q 3 g , and glueballs Q). A worldwide theoretical and 
experimental effort to search for dibaryon states with and 
without strangeness has lasted a long time already. 

Recently, an IJ P = 03 + deep bond A A resonance, 
called d* in 1989 fH, with a resonance mass M = 2.37 
GeV and a width r « 70 MeV had been experimentally 
confirmed by WASA-at COSY collaboration @-0]- This 
“inevitable” dibaryon has been extensively studied in our 
group fj, [§!, [l|. To provide more theoretical input for ex¬ 
perimental search, we included IJ P = 03 + A A channel 
coupling to the NN D- wave scattering phase shifts cal¬ 
culation and a resonance at 7? c . m . = 2.36 MeV had been 
shown in jToj | which called for the precise pn scattering 
measurement. Based on the fitting of our model to the 
NN and YN interaction data our research had also been 
extended to dibaryons with strangeness, several strange 
dibaryon resonances had been shown in specific baryon- 
baryon channels 0m, among which the most interest¬ 
ing one is the NQ state, a very narrow resonance in AH 
D-wave scattering. 

In 1987, Goldman et al. proposed that the S = —3, 
7=1/2, <7 = 2 dibaryon state NQ might be a narrow res¬ 
onance in a relativistic quark model [13j | . The quark de- 
localization color screening model (QDCSM) confirmed 
that it was a very narrow resonance 11], whereas the chi¬ 
ral quark model also claimed that TV12 might be bound 
state 0 | . Recently, the interest in the TV 12 dibaryon has 
been revived by lattice QCD calculations of HAL QCD 
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Collaboration [15]. They reported that the TV12 is in¬ 
deed a bound state at pion mass 875 MeV. To search 
for NQ dibaryon experimentally through the D —wave 
AH scattering process is complicated, as neither A nor 
H is a stable particle, H decays to An, and A decays to 
pn. Since N is stable, it is possible to observe NQ scat¬ 
tering experimentally. In order to provide more choices 
for experimental search, here we calculate the low-energy 
scattering phase shifts, scattering length, effective range 
and the binding energy of NQ system. If it is a sharp 
resonance, it could be observed by relativistic heavy-ion 
collisions data obtained at RHIC and LHC, or by the 
hadron beam experiments at J-PARC with their new de¬ 
veloped automatic scanning system |16j. 

Quantum chromodynamics (QCD) has been verified to 
be the fundamental theory of the strong interaction in the 
perturbative region. However, in the low energy region, 
it is hard to directly use QCD to study the complicated 
systems such as hadron-hadron interactions and multi¬ 
quark states due to the non-perturbative complication, 
although lattice QCD has made impressive progresses on 
nucleon-nucleon (NN) interactions and tetra- and penta- 
quark systems 0-0. Therefore, various QCD-inspired 
models have been developed to get physical insights of 
the multi-quark systems. 

To study the baryon-baryon interaction, the most com¬ 
mon approach is the chiral quark model (ChQM) |20l.l2l|. 
in which the constituent quarks interact with each other 
through colorless Goldstone bosons exchange in addition 
to the colorful one-gluon-exchange and confinement. To 
obtain the immediate-range attraction of TVTV interac¬ 
tion, the chiral partner cr meson-exchange has to be in¬ 
troduced. Although the a meson, as 7T7 t S'-wave reso¬ 
nance, had been observed by BES collaboration [22], the 
calculation of the nuclear force with correlated nn ex¬ 
change could not obtain enough attraction [23] as the 
phenomenological cr meson exchange did. 

An alternative approach to study baryon-baryon inter¬ 
action is the quark delocalization color screening model 
(QDCSM), which was developed in 1990s with the aim 
of explaining the similarities between nuclear and molec- 
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ular forces (13] • The model gives a good descrip¬ 
tion of NN and YN interactions and the properties of 
deuteron m. It is also employed to calculate the 
baryon-baryon scattering phase shifts in the framework 
of the resonating group method (RGM), and the dibaryon 
candidates are also studied with this model OKll. Re¬ 
cent studies also show that the NN intermediate-range 
attraction mechanism in the QDCSM, quark delocaliza¬ 
tion and color screening, is equivalent to the a meson- 
exchange in the chiral quark model, and the color screen¬ 
ing is an effective description of the hidden-color channels 
coupling [26l [27{| . 

So it is interesting to do a comparative study of the Nil 
system with these two quark models. This state might 
serve as a test of the flavor-dependent q — q interaction 
due to Goldstone-boson exchange and due to quark de- 
localization color screening, because there is no common 
flavor quark between N and f 1 and so no quark exchange 
between these two baryons. 

The structure of this paper is as follows. A brief in¬ 
troduction of constituent quark models used is given in 
section II. Section III devotes to the numerical results and 


discussions. The summary is shown in the last section. 


II. TWO QUARK MODELS 

In order to estimate the model dependence of dibaryon 
prediction, the two quark models: ChQM and QDCSM, 
are employed here to study the Nfl system. 


A. Chiral quark model 

The Salamanca model was chosen as the representa¬ 
tive of the chiral quark models, because the Salamanca 
group’s work covers the hadron spectra, nucleon-nucleon 
interaction, and multiquark states. In this model,the 
constituent quarks interact with each other through 
Goldstone boson exchange and one-gluon-exchange in ad¬ 
dition to the color confinement. The model details can 
be found in Ref. [21] . Here we only give the Hamiltonian: 
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Where a s is the quark-gluon coupling constant. In order 
to cover the wide energy scale from light, strange to heavy 
quark one introduces an effective scale-dependent quark- 
gluon coupling constant a s (fj,) HU , 


a s (n) 
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where g, is the reduced mass of the interacting quark-pair. 
The coupling constant g c h for chiral field is determined 
from the NNn coupling constant through 
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The other symbols in the above expressions have their 
usual meanings. 

For dibaryon with strangeness, Two versions of chi¬ 
ral quark model [H, [3Cj had been used. One is the so 
called extended chiral SU(2) quark model, in which a 
meson-exchange was used between any quark pair. An¬ 
other is the chiral SU(3) quark model, where full SU(3) 
scalar octet meson-exchange was used. In order to ana¬ 
lyze the effect of the scalar meson-exchange of the ChQM, 
in this work, we use three kinds of chiral quark models: 
(1) SU(2) ChQMl, where a meson-exchange is univer¬ 
sal, i.e., it exchanges between any quark pair; (2) SU(2) 
ChQM2, where a meson is restricted to exchange be¬ 
tween u and/or d quark pair only; (3) SU(3) ChQM, 
where the full SU(3) scalar octet meson-exchange is em¬ 
ployed. These scalar potentials have the same functional 
form as the one of SU(2) ChQM but a different SU(3) 
operator dependence [29|], that is, 
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B. Quark delocalization color screening model 

The Hamiltonian of QDCSM is almost the same as 
that of chiral quark model but with two modifications jl3, 
[25| : First, there is no er-meson exchange in QDCSM, and 
second, the screened color confinement is used between 
quark pairs reside in different baryon orbits. That is 

{ — a c A/ • A j (r^ + vo) if i,j in the same 

baryon orbit 

—a c A? ■ Xj ( 1-e ^ * J + v 0 ) otherwise 

( 12 ) 

The color screening constant fiij in Eq. (fl2l) is determined 
by fitting the deuteron properties, NN scattering phase 
shifts and N A, NT, scattering cross sections, n uu = 0.45, 
/ x us = 0.19 and /j, ss = 0.08, which satisfy the relation, 

k'us = k'uu * d’ss- 

The quark delocalization in QDCSM is realized by re¬ 
placing the left, right centered single Gaussian functions, 
the single particle orbital wave function in usual quark 
cluster model 
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= (<l> a (Si)+e<t> a (-Si))/N(e), 
M~Si,e) = (MSi) + eMSi))/N(e), 

N(e) = v / l + e 2 + 2ee-' s ?/ 4h2 . (14) 

The delocalization parameter e(Si) is determined by the 
dynamics of the quark system rather than adjusted pa¬ 
rameters. In this way, the system can choose its most 
favorable configuration through its own dynamics in a 
larger Hilbert space. 

The parameters of these models are given in Table |T| 
The calculated baryon masses in comparison with exper¬ 
imental values are shown in Table El 


TABLE I: The parameters of two models: m-n = 0.7fm _1 , 
mk = 2.51fm _1 , m v = 2.77fm _1 , m CT = 3.42fm _1 , m ao 
■ = rrif 0 — 4.97fm _1 , = 4.2fm _1 , A k = 5.2fm _1 , 

A v = 5.2fm _1 , A ct = 4.2fm _1 , A O0 = A K = A /o = 5.21m- 1 , 
5^/(470=0.54, 6>p=-15°. 



QDCSM 

SU(2)ChQMl SU(2)ChQM2 SU(3)ChQM 

6(fm) 

0.518 

0.518 

0.518 

0.518 

m u (MeV) 

313 

313 

313 

313 

md(MeV) 

313 

313 

313 

313 

m s (MeV) 

573 

573 

536 

573 

a c (MeV) 

58.03 

48.59 

48.59 

48.59 


0.45 

- 

- 

- 

Hus( fm -2 ) 

0.19 

- 

- 

- 

/r ss (fm -2 ) 

0.08 

- 

- 

- 

uo (MeV) 

-1.2883 

-1.2145 

-1.2145 

-0.961 

ao 

0.510 

0.510 

0.510 

0.583 

> 

o 

? 

1 

1.525 

1.525 

1.525 

1.616 

Mo (MeV) 

445.808 

445.808 

445.808 

422.430 


TABLE II: The masses of ground-state baryons(in MeV). 



N 

A 

A 

E 

E* 

H 

p—•* 

Q 

QDCSM 

939 

1232 

1124 

1238 

1360 

1374 

1496 

1642 

SU(2)ChQMl 

939 

1232 

1124 

1239 

1360 

1376 

1498 

1644 

SU(2)ChQM2 

939 

1232 

1137 

1245 

1376 

1375 

1506 

1620 

SU(3)ChQM 

939 

1232 

1123 

1267 

1344 

1398 

1475 

1625 

Expt 

939 

1232 

1116 

1193 

1385 

1318 

1533 

1672 


III. THE RESULTS AND DISCUSSIONS 

Here, we investigate the low-energy properties of 
the N-Q scattering with quantum numbers S, I, J = 
—3,1/2, 2 within various constituent quark models men¬ 
tioned above. The effects of channel-coupling are studied 
carefully, both color-singlet channels (the color symme¬ 
try of 3g-cluster is [111]) and hidden-color channels (the 
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color symmetry of 3<7-cluster is [21], which is called col¬ 
orful 3g-cluster) are included. The colorful 3g-cluster are 
listed in Table IIIII with color symmetry [c], spin symme¬ 
try [<r], flavor symmetry [/], isospin I, and strangeness 
S. The orbital symmetry is limited to be [3]. The labels 
of all 16 coupled channels of the Nil system are listed in 
Table |IV] 

To check whether or not there is a bound Nil state, 
a dynamic calculation based on the resonating group 
method (RGM) [31] has been done. Expanding the rel¬ 
ative motion wavefunction between two clusters in the 
RGM equation by Gaussian bases, the integro-differential 


equation of RGM reduces to algebraic equation-a gener¬ 
alized eigen-equation. The energy of the system is ob¬ 
tained by solving this generalized eigen-equation. In the 
calculation, the baryon-baryon separation is taken to be 
less than 6 fm (to keep the dimensions of matrix manage¬ 
ably small). The binding energies of Nil state in various 
quark models are listed in Table [V] where B sc stands for 
the binding energy of the single channel Nil, B§ cc means 
the binding energy with the five color-singlet channels 
included, and B ig cc refers to the binding energy with all 
the 16 channels coupling. 


TABLE III: The symmetries of colorful 3q— cluster. 



N' 

s' 

Tj/ 

A' 

E*' 

r-,*/ 

n' 

N” 

E” 

H 

A" 

a; 

[c] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

H 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[21] 

[3] 

[3] 

[3] 

[3] 

[21] 

[/] 

[21] 

[21] 

[21] 

[21] 

[3] 

[3] 

[3] 

[21] 

[21] 

[21] 

[21] 

[111] 

i 

1 

2 

1 

1 

2 

0 

1 

1 

2 

0 

1 

2 

1 

1 

2 

0 

0 

s 

0 

1 

2 

1 

1 

2 

3 

0 

1 

2 

1 

1 


TABLE IV: Channels of the Nil system. 


1 

2 

3 
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5 

6 

7 

8 
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EE* 

E*A 

Nil 


S”E*' 

E*'S” 

e” e” 

9 

10 

11 

12 

13 

14 

15 

16 

E”E 

N'il' 

E*'A” 

[i] 

> 

S” A” 

H”A' S 

E'E” 

E” A' 


The single channel calculation shows that the Nil is 
unbound (labelled as ’ub’ in Table [V]) in all quark models 
except the SU(2)ChQMl, in which a universal a meson- 
exchange is used. To show the contribution of each in¬ 
teraction term to the energy of the system, effective po¬ 
tentials between N and il are shown in Fig. 1, in which 
the contributions from various terms, the kinetic energy 
(Vufc), the confinement (V con ), the one-gluon-exchange 
{Voge), the one-boson-exchange (14, Vk, and V v ), and 
the scalar octet meson-exchange (14, 14 0 , 14, and Vf 0 ) 
are given. For the ChQM, the confinement, the one- 
gluon-exchange do not contribute to the effective poten¬ 
tial between N and il, because no quark exchange be¬ 
tween these two baryons, the pion and the cqO) meson 
do not contribute either because they do not exchange 
between u(d) and s quark. The contributions of other 
terms to the effective potential are shown in Fig. 1(a), 
Fig. 1(b) and Fig. 1(c), in which we can see that the 
attraction of the Nil system mainly comes from the a 
meson-exchange interaction. In SU(2)ChQMl, the at¬ 
traction from er meson-exchange is so large that it leads 
to a deep attractive potential between the N and il, which 


TABLE V: The binding energies B with channel-coupling. 



B sc (MeV) 

Bscc {MeV) 

Biecc {MeV) 

QDCSM 

ub 

-6.4 


SU(2)ChQMl 

-19.6 

-48.8 

-119.5 

SU(2)ChQM2 

ub 

ub 

-38.3 

SU(3)ChQM 

ub 

ub 

-13.7 


makes the Nil bound. While in SU(2)ChQM2, where the 
a meson is restricted to exchange between the u and d 
quarks only, there is no cr meson-exchange interaction be¬ 
tween N and il. So the total potential is repulsive, result¬ 
ing in unbound Nil. In SU(3)ChQM, even though the 
universal cr meson-exchange introduces large attraction, 
which is canceled by the repulsive potentials of n and fo 
meson-exchange, and also causes the Nil unbound. 

Things are different in QDCSM. In this model, quark 
delocalization and color screening work together to pro¬ 
vide short-range repulsion and intermediate-range attrac¬ 
tion. We illustrate this mechanism by showing contribu¬ 
tions of all interaction terms to the effective potential in 
Fig. 1(d), from which we see that the attraction of the 
Nil system mainly comes from the kinetic energy term 
due to quark delocalization, other terms provide repul¬ 
sive potentials, which reduce the total attraction of the 
Nil potential. The single channel approximation of the 
Nil in QDCSM is also unbound as shown in Table El 
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FIG. 1: The contributions to the effective potential from various terms of interactions. 


Then we consider the effects of channel-coupling. From 
Table lYl we can see that in SU(2)ChQMl, the binding en¬ 
ergy of the Nil increases to 48.8 MeV by including color- 
singlet channels coupling, and reaching the largest bind¬ 
ing energy of —119.5 MeV with the hidden-color channels 
coupling, which shows the channel-coupling has a huge 
influence on the Nil state. In the SU(2)ChQM2 and 
SU(3)ChQM, the effect of color-singlet channels coupling 
is not large enough to make the Nil bound, only addi¬ 
tional hidden-color channels coupling leads to the bound 
Nil state in both models. In QDCSM, since it contains 
hidden-color channels coupling effect already through the 
color screening EM, including the color-singlet chan¬ 
nels coupling only already makes the Nil state bound. 
All these results are consistent with our previous study of 
non-strange channels EE E3| , in which we found the QD¬ 
CSM with color-singlet channels copling only had similar 
results of the chiral quark models with both color-singlet 
and hidden-color channels coupling. 

If the Nil dibaryon is an A—wave bound state, the 
strong decays to S —wave octet-decuplet (£S*, S£*, and 
AS*) and decuplet-decuplet (£*S*) channels are prohib¬ 
ited kinematically because of its mass will be lower than 
the thresholds of these channels. It can only decay to 
the D —wave AS. We had done a scattering calculation 
by coupling the S'—wave Nil and the D —wave AS, and 
found the S —wave bound state Nil showed as a reso¬ 
nance in the D —wave AS scattering process El- Because 
of only the tensor interaction can couple the S —wave 
bound state to D —wave octet-octet baryon channel, so 
the energy of the bound state are pushed down only a lit¬ 
tle and the resonance width is generally small. We learnt 
from the Shanghai group of STAR collaboration that it 


is quite complicate to analyze the D —wave AS scattering 
data because both A and S are weak interaction unsta¬ 
ble.To provide more theoretical input, here we calculate 
the low-energy scattering phase shifts, scattering length 
and the effective range of Nil system. All results given 
below are calculated with 16 channels coupling in the chi¬ 
ral quark models and 5 color-singlet channels coupling in 
QDCSM. 

Firstly, we calculate the S'—wave Nil low-energy scat¬ 
tering phase shifts by using the well developed Kohn- 
Hulthen-Kato(KHK) variational method. The details 
can be found in Ref. [3l|. Fig. 2 illustrates the scat¬ 
tering phase shifts of the S —wave Nil. It is obvious that 
in all four quark models, the scattering phase shifts go to 
180 degrees at E c _ m . ~ 0 and rapidly decreases as E c . m , 
increases, which implies the existence of a bound state. 
The results are consistent with the the bound state calcu¬ 
lation shown before andthe lattice QCD calculation 0 - 

Secondly, by using these low-energy phase shifts, we 
extract the scattering length a o and the effective range 
r 0 of Nil scattering. Then the binding energy B' is ob¬ 
tained. The scattering length and effective range are cal¬ 
culated from the low-energy scattering phase shifts as: 

kcotS = -1- —rok 2 + 0(fc 4 ) (15) 

o-o 2 

The binding energy B' is calculated according to the re¬ 
lation: 


where /it is the reduced mass of the Nil; a is the wave 
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FIG. 2: The phase shifts of S —wave NQ dibaryon. 


number which can be obtained from the relation [321 ]: 

r-o = -(!-—) (17) 

a aao 

Please note that we use another method to calculate the 
binding energy, labelled as B'. The results are listed in 
Table ED 


TABLE VI: The scattering length ao, effective range r o, and 
binding energy B of the NQ dibaryon. 



a 0 (fm) 

r 0 (fm) 

B ( MeV) 

QDCSM 

2.8007 

0.5770 

-5.2 

SU(2)ChQMl 

0.8103 

0.3609 

-110.3 

SU(2)ChQM2 

1.3808 

0.6018 

-37.3 

SU(3)ChQM 

1.9870 

0.7064 

-13.7 


From Table eh we can see that in four quark models, 
the scattering length are all positive, which implies the 
existence of a bound state of NQ. The binding energies 
of the NQ dibaryon from the two methods (see B and 
B ) are coincident with each other. 


magnetic interaction. In this letter, we show that the N- 
Q system with quantum numbers SIJ = —3^2 their ef¬ 
fective interaction have very small contribution from the 
color-magnetic interaction due to its special quark con¬ 
tent. Because of this, the study of this state might teach 
us something about the mechanism of the intermediate- 
range attraction of the baryon-baryon interaction. In 
NN case we had shown that the phenomenological a 
meson exchange is equivalent to the quark delocaliza¬ 
tion and color screening and the latter is an effective de¬ 
scription of hidden color channel coupling. For the NQ 
system, the chiral quark model now only has the Gold- 
stone boson exchange, which will predict there is not a 
bound state if the er-meson is not universal exchange be¬ 
tween any quark pair. The physical observed a meson is 
an uu + dd system which should not exchange between 
u(d) and s quark. Only if one includes the hidden color 
channels coupling then the chiral quark model can ac¬ 
comodate a bound NQ state. On the other hand, the 
QDCSM predicts there will be a bound NQ state if the 
color singlet channels coupling has been taken into ac¬ 
count. The quark delocalization and color screening with 
the five color singlet channels coupling provides enough 
effective attraction to bound the NQ. Therefore if experi¬ 
ment confirms the existence of NQ dibaryon state, bound 
or appears as a D-wave resonance in A-S scattering, It 
will be a signal showing that the quark delocalization and 
color screening (an effective description of hidden color 
channels coupling) is really responsible for the intermedi¬ 
ate range attraction of baryon-baryon interaction. This 
mechanism is also preferred by the similarity between 
nuclear force and molecular force. 

Experimental confirmation of the NQ dibaryon will 
provide the second sample of six-quark system. It will 
provide another sample for the persue of low energy scale 
QCD property and we wish there will be more experimen¬ 
tal groups to be involved in the search of NQ dibaryon 
with quantum numbers SIJ P = —3^2 + . 
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